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Computational  fluid  dynamics  (CFD)  approaches  were  used  to 
compute  the  supersonic  and  h5^ersonic  flow  fields  and  aerod5mamic 
forces  and  moments  on  elliptical  projectiles.  Steady  state  numerical 
results  have  been  obtained  at  several  supersonic  Mach  niimbers 
between  2.5  and  4.0  and  several  angles  of  attack  from  0  to  12  degrees 
for  the  "jet-off"  conditions  with  the  use  of  Euler  and  Navier-Stokes 
flow  solvers.  In  addition,  numerical  computations  have  been 
performed  for  the  "jet-on"  conditions  to  study  the  interaction  of  a 
helixim  jet  with  a  free  stream  Mach  4.0  flow.  In  general,  very  good 
agreement  of  the  computed  aerodynamic  coefficients  with  the 
experimental  data  was  achieved  at  all  speeds  and  angles  of  attack 
investigated  for  jet-off  conditions.  A  small  discrepancy  exists  in  the 
comparisons  for  the  axial  force.  CFD  results  for  the  jet-on  case  show 
the  qualitative  features  and  strong  flow  interaction  between  the  jet 
and  the  free  stream  flow.  The  results  show  the  predictive  capabilities 
of  CFD  techniques  for  supersonic  flow  over  elliptical  projectiles. 
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NUMERICAL  COMPUTATIONS  OF  SUPERSONIC  FLOW 
OVER  ELLIPTICAL  PROJECTILES 


1.  Introduction 


The  advancement  of  computational  fluid  d)mamics  (CFD)  has  had  a  major 
impact  on  projectile  design  and  development  [1,2^4].  Improved  computer 
technology  and  state-of-&e-art  numerical  procedtires  enable  solutions  to 
complex,  three-dimensional  (3-D)  problems  associated  with  projectile  and  missile 
aerodynamics.  In  general,  these  techniques  produce  accurate  and  reliable 
numerical  results  for  projectiles  and  missiles  at  small  angles  of  attack.  Modem 
projectiles  and  missiles  are  expected  to  experience  moderate  to  large  angles  of 
attack  during  flight.  Of  partictdar  interest  is  the  accurate  determination  of 
supersonic  and  hypersonic  flow  over  elliptical  projectiles  at  moderate  angles  of 
attack.  The  flow  field  for  such  projectiles  with  non-axisymmetric  cross  sections  is 
complex,  especially  in  the  presence  of  jets  used  to  maneuver  these  projectiles. 
The  work  presented  in  this  report  was  initiated  ais  part  of  The  Technical 
Cooperation  Program  (TTCP)  effort  with  participants  feom  Canada,  the  United 
Kingdom,  and  the  United  States  and  was  aimed  at  assessing  the  capabilities  of 
the  both  Euler  and  Navier-Stokes  solvers  currently  available  to  research  scientists 
for  supersonic  and  hypersonic  flow  over  elliptical  projectiles  for  both  "jet-off" 
and  "jet-on"  conditions  [5].  The  TTCP  research  effort  has  also  focused  on  the 
wind  tunnel  testing  as  well  as  free  flight  testing  of  these  projectiles.  Different 
aspects  of  computational  techniques,  such  as  grid  generation,  algorithms, 
turbtilence  modeling  and  flow  field  visualization,  have  been  addressed  by  the 
group. 

Inviscid  solutions  were  obtained  for  H-series  projectiles  with  the  use  of  the  Euler 
Zeus  graphical  user  interface  (GUI)  code  [6,7].  Computations  of  the  H3  projectile 
and  its  variations  (with  and  without  flares  and  strakes)  were  performed  with  the 
Euler  Zeus  code  at  M  =  8.2  and  several  angles  of  attack  between  0®  and  15®. 
Comparison  was  made  of  stabilization  by  strakes  and  flare  for  a  different  H-series 
projectile  [8].  The  computations  indicated  that  altiiough  the  flare  gave  increased 
drag,  it  was  still  substantially  more  effective  in  providing  stability  than  strakes 
alone.  The  present  research  focuses  on  the  application  of  CFD  techniques  for 
accurate  numerical  prediction  of  supersonic  flow  over  the  elliptical  H3P78 
projectile.  Calculations  for  the  H3P78  projectile  were  performed  with  the  Zeus 
Euler  code  and  two  Navier-Stokes  flow  solvers:  the  zonal  Navier-Stokes  flow 
(solver)  (ZNSFLOW)  [9]  and  CFD++  [10,11],  at  several  supersonic  Mach  numbers 
between  2.5  and  4.0  and  several  angles  of  attack  from  0®  to  12®  for  flie  jet-ofl 
conditions.  In  addition,  numerical  computations  for  the  H3P78  projectile  have 
been  performed  for  the  jet-on  conditions  with  the  CFD++  code  to  study  Ihe 
interaction  of  a  helium  jet  with  a  free  stream  M  =  4.0  flow. 


1 


A  description  of  the  computational  techniques  is  presented,  followed  by  a 
description  of  the  applications  of  these  techniques  to  the  H3P78  projectile. 
Results  for  this  configuration  are  shown  at  various  supersonic  speeds  and  angles 
of  attack.  Computed  data  have  been  compared  witii  experimental  data  provided 
by  the  Defence  Evaluation  and  Research  Agency  (DERA),  United  Kingdom  (UK) 
[51  and  obtained  at  the  Defence  Research  Establishment,  Valcartier  (DREV), 
Canada  [12]. 


2.  Solution  Technique 


2.1  Euler  Zeus  GUI  Code 

Development  of  the  Zeus  GUI  code  [6,7]  at  the  U.S.  Army  Research  Laboratory 
(ARL)  was  undertaken  at  tiie  request  of  die  Aviation  and  Missile  Research 
Development  and  Engineering  Center  (AMRDEC)  and  was  fimded  by  ARL  as 
part  of  a  Technology  Program  Annex  (TP A)  agreement.  A  picture  of  the  Zeus 
GUI  code  in  use  is  shown  in  Figure  1.  As  a  starting  point  for  this  effort,  the  zonal 
Euler  solver  or  Zeus  code  [13,14]  was  chosen  by  AMRDEC  to  be  incorporated 
into  a  GUI  environment.  The  original  Zeus  code  was  developed  by  the  Navy  and 
is  widely  used  within  the  international  missile  design  community.  Zeus  is  a 
zonal  Elder  CFD  solver,  which  employs  a  second  order  Godimov  scheme  to 
integrate  the  Euler  equations  and  march  the  solution  longitudinally  along  the 
body.  The  solver  algorithm  in  the  ARL  Zeus  GUI  is  identical  to  the  one 
implemented  in  the  original  Zeus  code.  One  can  obtain  additional  details  about 
the  integration  scheme  by  referring  to  the  work  of  Wardlaw  [13,14]. 

In  terms  of  the  application  of  Zeus,  there  are  some  restrictions  of  its  use.  First, 
Zeus  can  only  be  applied  to  cases  in  which  the  entire  flow  field  is  supersonic. 
Zeus  employs  the  Euler  equations  and  tiierefore  is  most  suitable  for  cases  in 
which  the  boundary  layer  is  assumed  to  be  thin  and  its  effects  can  be  neglected. 
The  computational  mesh  should  be  free  of  blunt  discontinuities.  Auxiliary 
programs  are  available  to  generate  Zeus-compatible  solution  planes  for  blunt 
noses  when  the  computational  mesh  is  not  compatible  with  Zeus  topology  and 
when  subsonic  flow  may  occur.  Zeus  was  written  to  support  a  zonal  grid 
topology.  The  zonal  topology  allows  fins  or  wings  to  be  modeled.  However,  the 
leading  and  trailing  edges  should  be  fairly  sharp. 


2 


Figure  1.  Zeus  GUI. 

2.2  ZNSFLOW  Solver 

The  ZNSFLOW  solver  is  a  product  of  a  common  high  performance  computing 
software  support  initiative  (CHSSI)  project  [9].  It  is  a  descendant  of  F3D  [16],  a 
code  used  successfully  for  many  years  on  Cray  vector  processors  such  as  the  C90. 
Under  CHSSI,  it  was  rewritten  to  provide  scalable  performance  on  a  number  of 
computer  architectures.  Programming  enhancements  include  the  use  of  dynamic 
memory  allocation  and  highly  optimized  cache  management.  ZNSFLOW 
features  a  GUI  to  facilitate  problem  setup.  It  has  been  used  extensively  in  the 
computation  of  flow  field  calculations  for  projectile  and  missile  programs  of 
interest  to  the  U.S.  Army  [1,2,3,4,9,15].  The  flow  solver  includes  the  chimera 
[17,18,19]  overset  discretization  technique  for  CFD  modeling  of  complex 
configurations.  By  using  the  chimera  technique,  one  can  greatly  simplify  the  grid 
topology  and  grid  generation  for  very  complex  systems.  One  of  the 
disadvantages  in  using  the  chimera  technique  has  been  the  increased  complexity 
and  corresponding  confusion  in  the  application  of  a  turbulence  model.  A 
chimera  model  can  be  composed  of  multiple  zones,  with  each  zone  possibly 
having  a  unique  grid  topology.  Most  turbulence  models  have  specific  direction-, 
orientation-,  or  distance-related  requirements  for  correct  application.  For  a 
complex  chimera  model,  appl5dng  a  turbulence  model  can  be  a  very  complex 
process,  especially  an  algebraic  model  such  as  Baldwin-Lomax  [20]  model.  We 
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have  addressed  this  problem  in  ZNSFLOW  by  installing  point-wise  one-  and 
two-equation  turbulence  models  [21]  that  are  not  orientation  specific.  This 
greatly  simplifies  the  setup  of  the  turbulence  model.  Wall  location  information  is 
supplied  when  the  wall  bormdary  conditions  are  set  by  the  user. 

2.3  Governing  Equations 

The  complete  set  of  3-D,  time-dependent,  generalized  geometry,  Reynolds- 
averaged,  ihin  layer  Navier-Stokes  equations  is  solved  numerically  to  obtain  a 
solution  to  this  problem  and  can  be  written  in  general  spatial  coordinates  4,  r], 
and  ^  as  follows  [22]: 

d^q  +  d^F  +  d^G  +  d^H  =  RQ-^d^S,  (1) 

in  which 

^  =  4  (x,  y,  z,  t)  -  longitudinal  coordinate, 
ri  =  Ti  (x,  y,  z,  t)  -  circumferential  coordinate, 

=  (^  (x,  y,  z,  t)  -  nearly  normal  coordinate,  and 
T  =  t  -  time. 


In  Equation  1,  q  contains  the^  dependent  variables  (density,  three  velocity 
components,  and  energy),  andF,  G,  and  H  are  flux  vectors.  The  thin  layer 
approximation  is  used  here,  and  the  viscous  terms  involving  velocity  gradients 
in  both  the  longitudinal  and  circumferential  directions  are  neglected.  The  viscous 
terms  are  retained  in  the  normal  direction,  ^  and  are  collected  into  the  vector,  S . 
In  the  wake  or  the  base  region,  similar  viscous  terms  [1]  are  also  added  in  the 
streamwise  direction,  An  implicit,  approximately  factored  scheme  is  used  to 
solve  these  equations.  For  the  computation  of  turbulent  flows,  the  turbulent 
contributions  are  supplied  through  an  algebraic  eddy  viscosity  turbulence  model 
developed  by  Baldwin  and  Lomax  [20]  or  a  point-wise  turbulence  model  [21]. 

2.4  Numerical  Algorithm 

The  implicit,  approximately  factored  scheme  for  the  Ihin  layer  Navier-Stokes 
equations  with  central  differencing  in  the  t|  and  ^  directions  and  ^^upwinding”  in 
^  is  written  in  the  following  form  [16]: 

[  I  +  i.hS^^  iry+hhS^C"-i,hRQ-^  Si  ] 

x[  I  +  ithSy  (A-r+i,hS,B"-i,D,  1^  ]  AQ” 

=  4At  {  s!  [  (F^r-P:  ]  +  ^/  [(F-)"-F;]  +  ^,  (2) 

+  (H^-HJ-Re-'S,(S".S„)  }-4W-eJ  , 
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in  which  h  =  At  or  (At)/ 2  and  the  free  stream  base  solution  is  used.  The  free 
stream  fluxes  are  subtracted  from  the  governing  equation  to  reduce  the 
possibility  of  error  from  the  free  stream  solution,  which  corrupts  the  converged 
solution.  Here,  6  is  typically  a  three-point  second  order  accurate  central 
difference  operator,  <5  is  a  midpoint  operator  used  with  the  viscous  terms,  and 
the  operators  <5/  arid  S/  are  backward  and  forward  three-pomt  difference 
operators.  The  flux  F  has  been  eigensplit,  and  the  matrices  A,  B,C,  and  M 
resiflt  from  local  linearization  of  the  fluxes  about  the  previous  time  level.  Here, } 
denotes  the  Jacobian  of  the  coordinate  transformation.  Dissipation  operators  De 
and  D,  are  used  in  the  central  space-differencing  directions.  The  smoothing  terms 
used  in  the  present  study  are  of  the  form 


and 


DJ,  =  (A0^-* 


£2  Sp{E) P 5  +€/^^  5  ^ 

l  +  P 


I, 


A L  =  (A 0  J-'  [ f 2 S piB) pS +2.5 £,5p{B)5  ]  |,  J . 


in  which 


fP\ 


\ii  +  sn  p 


(3) 


(4) 


and  in  which  p(B)  is  the  true  spectral  radius  of  B.  The  idea  here  is  that  the  fourth 
difference  will  be  tuned  down  near  shocks  (e.g.,  as  p  gets  large,  the  weight  on  the 
fourth  difference  drops  down  while  the  second  difference  tunes  up). 

2.5  CFD++  Flow  Solver 

The  basic  numerical  framework  in  which  the  proposed  scheme  is  implemented  is 
termed  tiie  unified  grid,  unified  physics  and  unified  computing  framework. 
These  have  been  implemented  in  a  software  suite  called  CFD++  [10,11]  and  the 
user  is  referred  to  these  references  for  details  of  the  basic  numerical  framework. 
Here,  only  a  brief  S5mopsis  of  this  frcunework  and  methodology  is  given. 


The  3-D,  time-dependent  Reynolds-averaged  Navier-Stokes  (RANS)  equations 
are  solved  by  the  finite  volume  method: 

—  fWdF  +  -  G]-  =  \lAdV  (5) 

dt  y  y 

in  which  W  is  tihe  vector  of  conservative  variables,  F  and  G  are  the  inviscid  and 
viscous  flux  vectors,  respectively,  H  is  the  vector  of  source  terms,  V  is  the  cell 
volume,  and  A  is  the  surface  area  of  the  cell  face. 
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The  nximerical  framework  of  CFD++  is  based  on  the  following  general  elements: 

1.  Unsteady  compressible  and  incompressible  Navier-Stokes  equations 
with  turbulence  modeling  [unified  physics]; 

2.  Unification  of  Cartesian,  structured  curvilinear,  and  unstructured 
grids,  including  hybrids  [unified  grid]; 

3.  Unification  of  treatment  of  various  cell  shapes  including  hexahedral, 
tetrahedral,  and  triangular  prism  cells  (3-d),  quadrilateral  and  triangular  cells  (2- 
d),  and  linear  elements  (1-d)  [tmified  grid]; 

4.  Treatment  of  multi-block  patched  aligned  (nodally  connected), 
patched  nonaligned  and  overset  grids  [unified  grid]; 

5.  Total  variation  diminishing  discretization  based  on  a  new  multi¬ 
dimensional  interpolation  framework; 

6.  Riemann  solvers  to  provide  proper  signal  propagation  physics, 
including  versions  for  preconditioned  forms  of  the  governing  equations; 

7.  Consistent  and  accurate  discretization  of  viscous  terms  via  the  same 
multi-dimensional  polynomial  framework; 

8.  Point-wise  turbulence  models  that  do  not  require  knowledge  of 
distance  to  walls; 

9.  Versatile  boimdary  condition  implementation  includes  a  rich  variety 
of  integrated  boxmdary  condition  types  for  the  various  sets  of  equations;  and 

10.  Implementation  on  massively  parallel  computers  based  on  the 
distributed  memory  message-passing  model  via  native  message-passing  libraries 
or  message-passing  interface,  parallel  virtual  machines,  etc.  [tmified  computing]. 

The  code  has  assembled  several  ideas  about  convergence  acceleration  to  yield  a 
fast  methodology  for  all  flow  regimes.  The  approach  can  be  labeled  as  a 
"preconditioned  implicit  relaxation"  scheme.  It  combines  three  bcisic  ideas: 
implicit  local  time  stepping,  relaxation,  and  preconditioning.  Preconditioning  the 
equations  ideally  equalizes  the  eigenvalues  of  the  inviscid  flux  Jacobians  and 
removes  the  stiffness  arising  from  large  discrepancies  between  the  flow  and 
sound  velocities  at  low  speeds.  The  use  of  an  implicit  scheme  circumvents  die 
stringent  stability  limits  suffered  by  their  explicit  coimterparts,  and  successive 
relaxation  allows  cell  revision  as  information  becomes  available  and  thus  aids 
convergence. 

The  code  has  recently  added  the  ability  to  handle  multi-block  meshes  widi 
various  types  of  inter-block  connectivities.  Multi-dimensional  interpolation  more 
accurately  represents  local  behavior  of  flow-dependent  variables.  While  the 
formal  order  of  accuracy  does  not  need  to  be  any  higher,  this  approach  leads  to 


6 


practically  higher  accuracy  on  relatively  coarse  meshes.  The  multi-dimensional 
interpolation  framework  helps  us  to  easily  handle  inter-block  connectivities  also. 
Second  order  discretization  was  used  for  the  flow  variables  and  the  turbulent 
viscosity  equation.  The  turbulence  closure  has  been  based  on  topology- 
parameter-free  formulations.  These  models  are  ideally  suited  to  unstructured 
bookkeeping  and  massively  parallel  processing  because  of  their  independence 
from  constraints  related  to  the  placement  of  boundaries  and/or  zonal  interfaces. 
Recent  contributions  to  these  models  include  (a)  improved  behavior  of  the 
dissipation  rate  transport  equation  by  explicit  sensitization  to  non-equilibrium 
flow  regions,  and  (b)  enhanced  near- wall  characteristics  and  elimination  of  ad  hoc 
formulations  through  the  introduction  of  time-scale  realization. 


3.  Chimera  Composite  Grid  Scheme 


The  chimera  overset  grid  technique  greatly  adds  to  the  number  of  applications  to 
which  the  ZNSFLOW  and  CFD++  solvers  can  be  applied.  The  chimera  overset 
grid  technique,  which  is  ideally  suited  to  complex  configurations  and  multi-body 
problems  [1,2,3,4,17,18,19],  involves  generating  independent  grids  about  each 
body  or  component  and  then  oversetting  tiiem  onto  a  base  grid  to  form  tihe 
complete  model.  An  advantage  of  the  overset  grid  technique  is  that  it  allows 
computational  grids  to  be  obtained  for  each  body  component  separately  and  thus 
makes  the  grid  generation  process  easier.  Because  each  component  grid  is 
generated  independently,  portions  of  one  grid  may  lie  within  a  solid  boundary 
contained  within  another  grid.  Such  points  lie  outside  the  computational  domain 
and  are  excluded  from  the  solution  process.  Equation  2  has  been  modified  for 
chimera  overset  grids  by  the  introduction  of  the  flag  ibto  achieve  just  that.  This  ib 
array  accommodates  the  possibility  of  having  arbitrary  holes  in  the  grid.  The  ib 
array  is  defined  so  that  i  b  =  1  at  normal  grid  points  and  i  b  =  0  at  hole  points. 
Thus,  when  ib  =  1,  Equation  2  becomes  the  standard  scheme.  The  set  of  grid 
points  that  forms  tiie  border  between  the  hole  points  and  the  normal  field  points 
is  called  inter-grid  boimdary  points.  We  revise  these  points  by  interpolating  the 
solution  from  the  overset  grid  ttiat  created  the  hole.  Values  of  die  ib  array  and  the 
interpolation  coefficients  needed  for  this  revision  are  provided  by  a  separate 
algorithm  [17J. 


4.  Model  Geometry  and  Computational  Grid 


4.1  Zeus  Computational  Mesh 

One  of  the  more  convenient  features  of  the  Zeus  GUI  is  the  automatic  grid 
generation.  The  user  merely  specifies  the  missile  surface  geometry,  and  die  flow 
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field  mesh  is  built  plane  by  plane  as  the  solution  is  computed.  The  current  set  of 
GUI-directed  geometry  input  is  capable  of  describing  axisymmetric  finned 
missiles.  Since  H-series  projectiles  are  not  axis5munetric,  the  projectile  surface 
was  generated  by  an  auxiliary  program  and  written  as  a  set  of  discrete  points  to  a 
file.  The  projectile  surface  was  read  into  Zeus  GUI  via  the  user-defined  surface 
option.  The  Zeus  GUI  also  provides  controls  for  mesh  clustering.  The  mesh  was 
clustered  in  the  direction  normal  to  the  body  at  the  surface.  Pitch-plane 
symmetry  was  used  for  most  computations.  When  pitch-plane  symmetry  was 
employed,  72  points  were  used  in  the  circumferentid  direction.  The  outer  grid 
bormdary  of  the  automaticadly  generated  grid  is  shock  fitted.  Therefore,  the 
distance  from  die  projectile  body  surface  to  the  outer  grid  body  grows  as  the 
solution  is  marched  from  nose  to  tail.  In  most  cases,  36  points  were  used  in  the 
normal  direction  for  the  first  15%  of  body  length;  then  the  solution  was  stopped. 
The  flow  field  solution  for  the  remaining  body  length  was  usually  computed 
with  72  points  in  the  direction  normal  to  the  body,  but  for  some  high  angle-of- 
attack  cases,  120  points  were  used.  To  increase  the  number  of  points  in  the 
normal  direction,  an  auxiliary  program  was  used  to  interpolate  a  solution  plane 
with  36  normal  points  to  one  with  72  normal  points.  The  auxiliary  program  is 
accessible  through  the  GUI  and  uses  the  GUI  parameters,  such  as  mesh  cluster 
controls,  to  easily  refine  the  computational  mesh.  At  this  mesh  resolution,  Zeus 
usually  required  approximately  5,000  axial  planes  to  obtain  the  solution  for  the 
entire  body  length.  A  solution  could  be  obtained  with  approximately  20  minutes 
of  computer  time  on  a  Silicon  Graphics  workstation. 

4.2  Navier-Stokes  Computational  Meshes 

For  the  jet-off  cases  that  use  the  Navier-Stokes  flow  solvers,  two  multi-block 
grids  were  developed.  These  structured  multi-block  grids  have  one-to-one 
overlaps  at  the  zonal  boundaries.  Both  grids  are  two  zone  H-gnds.  The  H3P78 
projectile  is  shown  in  Figure  2.  The  initial  grid  (GRID-1)  used  for  these 
computations  was  a  two-zone  grid  (see  Figure  3)  that  consisted  of  1.8  milhon 
grid  points.  Zone  1  grid  along  the  projectile  body,  has  251  longitudinal  points,  59 
normal  points,  and  91  circumferential  points,  with  an  FI-type  grid  at  the  nose  of 
the  projectile.  The  base  grid.  Zone  2,  is  an  H-gnd  consisting  of  50  longitudinal 
points,  113  normal  points,  and  91  circumferential  points.  The  minimum  spacing 
at  the  wall  is  6.0E-05  mm  (lE-06  calibers).  An  expanded  view  of  the  base  region 
grid  is  shown  in  Figure  4. 


Consisting  of  one  million  points,  GRID-2  contains  fewer  total  points  but  has  a 
larger  number  of  points  (121  points,  see  Figure  5)  in  the  circumferential  direction 
in  both  zones.  Zone  1  has  101  longitudinal  points  and  45  normal  pointe,  while 
Zone  2,  in  the  wake  region,  has  40  longitudinal  points  and  99  normal  points.  The 
minimum  spacing  at  the  projectile  wall  is  2.0E-03  mm  (3.0E-05  calibers)  for  both 
grids. 
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82.8  mm 


*  Reference  area  for  forces  is  based  on 
projected  base  area  for  3/4  power  law 
surface  extended  to  the  base.  This 
would  be  60  mm  in  width.  .. 


Figure  2.  H3P78  Projectile. 
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Figure  4.  Expanded  View  of  the  Grid  in  the  Base  Region,  Jet  Off. 


Figure  5.  Circumferential  Cross  Section  of  the  Grid  at  the  Base. 
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For  the  jet-on  cases,  similar  multi-block  grids  for  the  body  were  generated.  The 
full  grid  for  the  jet-on  case  is  shown  in  Figure  6,  while  an  expanded  view  of  tire 
base  region  grid  is  shown  in  Figure  7.  Here,  for  ease  of  grid  generation,  tire 
chimera  overset  gridding  technique  (described  earlier)  was  used  to  model  the  jet. 
The  projectile  body  grid  was  generated  first.  The  chimera  technique  allows  the  jet 
to  be  gridded  separately.  Two  additional  grid  zones  were  created  to  model  the 
jet.  The  dimensions  for  each  zone  are  21  x  20  x  40  and  10  x  10  x  40,  in  the 
longitudinal,  circumferential,  and  normal  directions.  These  two  zones  were  then 
added  to  the  projectile  grid  for  the  jet-on  calculations,  adding  about  20,000  points 
to  each  grid.  One  of  the  jet  grids  is  a  cylindrical  grid  that  covers  the  actual  jet  and 
extends  beyond  (see  Figure  8).  The  second  one  is  a  rectangular  grid  placed 
within  die  actual  jet  and  overset  onto  the  other  grid  primarily  to  avoid  die  grid 
singularity  along  the  centerline  of  the  jet.  An  expanded  view  of  the  jet  grids 
projected  onto  the  projectile  surface  is  shown  in  Figure  8.  The  jet  grids  and  the 
body  grids  are  aU  overset  to  form  the  complete  mesh  system.  The  chimera 
procedure  results  in  hole  boundaries  (not  shown  here)  in  the  body  grid,  which 
are  attributable  to  the  jet,  and  transfers  information  between  the  jet  grids  along 
these  hole  boundaries.  The  outer  boundaries  of  the  jet  grids  also  receive 
information  interpolated  from  the  body  mesh. 


Figure  6.  FuU  Grid  for  Viscous  Computations,  Jet  On. 


A  fine  mesh  (GRID-3)  was  developed  for  the  jet-on  calculations  (see  Figure  6). 
With  GRID-2  as  a  starting  point,  the  number  of  longitudinal  points  was  increased 
to  140,  with  the  points  more  strategically  placed  along  the  projectile  flare  in  the 
vicinity  of  the  jet.  The  other  grid  zones  were  unchanged  and  were  the  same  as 
previously  described  for  GRID-2.  The  total  number  of  points  for  GRID-3  is  1.3 
million.  Again,  the  minimum  spacing  at  the  wall  is  2.0E-03  mm. 
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Figure  7.  Expanded  View  of  the  Grid  in  the  Base  Region,  Jet  On. 
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5.  Results 


Steady  state  numerical  computations  using  botii  Euler  and  viscous  Navier-Stokes 
methods  were  performed  to  predict  the  flow  field  and  aerodynamic  coefficients 
on  the  H3P78  elliptical  projectile  for  both  jet-off  and  jet-on  conditions.  Three- 
dimensional  calculations  were  performed  at  supersonic  speeds  between  Mach  2.5 
and  4.0  and  several  angles  of  attack  between  0°  and  12°.  The  projectile  geometry 
and  a  set  of  experimental  wind  tunnel  data  for  validation  of  the  computations 
were  supphed  by  DERA,  UK,  and  DREV,  Canada. 

5.1  Zeus  Euler  Code  Results  (jet  off) 

Presented  are  results  of  computations  for  the  H3P78  configuration.  The  Zeus  GUI 
code  was  used  to  compute  the  flow  field  of  the  H3P78  at  Mach  numbers  2, 3,  and 
4.  For  each  Mach  number,  the  flow  field  was  computed  at  multiple  angles  of 
attack.  Figures  9, 10,  and  11  show  pressure  contours  of  the  H3P78  body  and  flow 
field  at  Mach  4.  For  the  most  part,  Zeus  captures  the  flow  field  features  fairly 
weU.  However,  Zeus  was  shown  to  be  more  accurate  in  predicting  force  and 
moments  at  lower  angles  of  attack  than  at  higher  angles  of  attack  for  H-series 
projectiles.  This  may  indicate  that  the  lee  side  separation  at  higher  angles  of 
attack  is  producing  flow  field  features  that  Zeus  is  not  predicting  accurately. 
Also  note  that  Zeus  can  only  perform  the  forebody  flow  field  computation. 
Therefore,  force  and  moments  do  not  include  the  contribution  of  the  base  flow 
field. 


Figure  9.  Pressure  Contours,  M  =  4,  a  =  4°. 
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Figure  10.  Pressure  Contours,  M  =  4,  a  =  8°. 


Figure  11.  Pressure  Contours,  M  =  4,  a  =  12°. 


Force  and  moments  computed  from  the  H3P78  flow  fields  are  shown  in 
Figures  12  through  15.  In  feese  cases,  the  moment  coefficients  were  calculated 
about  the  nose  of  the  projectile.  Perhaps  the  most  useful  information  that  can  be 
determined  from  a  review  of  the  plots  is  the  information  in  Figure  15  that  shows 
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how  the  forces  are  generated  with  body  length.  The  plot  indicates  the  extent  to 
which  the  flare  of  the  H3P78  contributes  to  the  overall  stability  of  the  projectile. 
In  general,  the  flare  contributes  greatly  to  the  stability,  but  Figure  15  also 
indicates  the  substantial  drag  penalty  attributable  to  the  flare.  This  is  indicated  in 
Figure  15  by  the  increase  in  the  slope  of  tire  axial  force  cxxrve.  In  another  study 
performed  with  a  Zeus  GUI,  a  comparison  was  made  of  stabilization  by  strokes 
and  flare  for  a  different  H-series  projectile  [8].  The  computations  indicated  that 
although  the  flare  gave  increased  drag,  it  was  still  substantially  more  effective  in 
providing  stability  than  strakes  alone.  Overall,  Zeus  is  a  very  good  design  tool  for 
analyzing  this  type  of  trade-off. 


Figure  12.  Axial  Force  Versus  Angle  of  Attack. 


Figure  13.  Normal  Force  Versus  Angle  of  Attack. 
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Figure  14.  Pitching  Moment  Versus  Angle  of  Attack. 
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Figure  15.  Force/ Moment  Versus  Body  Length. 
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5.2  Navier-Stokes  Results  (jet  off) 

Three-dimensional  numerical  computations  have  been  performed  for  the  H3P78 
projectile  with  both  ZNSFLOW  and  CFD++  codes  at  several  supersonic  Mach 
numbers  between  2.5  and  4.0  and  several  angles  of  attack  from  0°  to  12°  for  the 
jet-off  conditions.  In  addition,  numerical  computations  for  the  H3P78  projectile 
have  been  performed  for  the  jet-on  conditions  to  study  tiie  interaction  of  a 
helium  jet  with  a  free  stream  M  =  4.0  flow.  All  these  simulations  employ 
symmetry  in  the  pitch  plane,  so  only  half  of  the  physical  domain  is  modeled. 
Turbulence  closure  was  provided  via  either  an  algebraic  model  or  a  point-wise 
two-equation  model.  Computed  results  obtained  for  the  jet-off  conditions  are 
presented  next. 

Figures  16  and  17,  respectively,  show  computed  pressure  contours  at  M  =  2.5  and 
4.0  for  different  angles  of  attack.  Here,  blue  represents  low  pressure  and  red 
represents  high  pressure  regions.  Figure  16  shows  the  computed  pressure 
contours  at  Mach  2.5  and  angles  of  attack  4,  8,  and  12  degrees  (from  top  to 
bottom).  It  shows  the  higher  pressures  on  the  wind  side  and  lower  pressures  on 
tire  lee  side.  The  lower  pressures  on  the  lee  side  can  be  seen  more  clearly, 
especially  at  12  degrees  angle  of  attack.  Also  seen  in  this  figure  are  the  shock 
waves  emanating  from  the  nose  and  the  body-flare  jtmction.  Similar  flow  field 
features  can  be  seen  at  Mach  4.  As  seen  in  Figure  17,  the  shock  waves  are 
stronger  for  M  =  4  both  at  the  nose  and  the  flare  jrmction.  As  the  angle  of  attack 
increases  from  4  to  12  degrees  (from  top  to  bottom),  the  computed  pressure 
contours  show  an  increase  in  the  low  pressure  region  on  tiie  lee  side.  Figure  18 
shows  the  computed  particle  traces  in  the  base  region  at  M  =  4.0  and  a  =  12 
degrees.  These  particle  traces  show  the  complex  vortical  flow  pattern  in  the  near 
wake  of  the  projectile. 

Figures  19,  20,  and  21  show  the  comparison  of  the  computed  aerodynamic  force 
and  moment  coefficients  with  the  data  at  M  =  2.5.  We  obtained  these  computed 
force  and  moment  coefficients  by  integrating  Ihe  computed  surface  pressures. 
Figure  19  shows  the  comparison  of  the  computed  normal  force  coefficient  with 
the  data.  The  computed  results  shown  here  were  obtained  witii  both  ZNSFLOW 
and  CFD++  flow  solvers.  ZNSFLOW  results  were  obtained  with  two  different 
computational  grids  (grid  1  consisting  of  1.8  million  and  grid  2  consisting  of 
1.0  million  grid  points)  and  show  tiiat  ttie  computed  results  are  grid 
independent.  The  computed  normal  force  obtained  by  botii  ZNSFLOW  and 
CFD++  codes  matches  fairly  well  with  the  data  except  at  12  degrees  angle  of 
attack  where  a  small  discrepancy  exists.  The  computed  results  with  the 
ZNSFLOW  code  were  obtained  with  an  algebraic  turbulence  model,  whereas  the 
CFD++  results  were  obtained  with  a  point-wise  two-equation  turbulence  model. 
Figure  20  shows  the  comparison  of  the  computed  pitching  moment  coefficient 
with  the  experimental  data.  Again,  the  nose  of  the  projectile  is  used  as  Ihe 
reference  point  for  the  moment  coefficient.  In  general,  the  computed  pitching 
moment  coefficient  obtained  by  both  Navier-Stokes  solvers  is  in  good  agreement 
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with  the  data,  regardless  of  the  turbulence  modeling  used  in  the  computations. 
The  data  and  the  computed  results  show  that  die  pitching  moment  coefficient 
decreases  as  the  angle  of  attack  increases.  A  comparison  of  the  computed  axial 
force  coefficient  witii  the  data  is  shown  in  Figure  21.  The  data  and  the  computed 
resxilts  show  that  the  drag  is  almost  constant  as  a  fimction  of  angle  of  attack. 
Again,  the  computed  axial  force  coefficients  are  in  good  agreement  and  are 
within  10%  of  the  measured  data.  The  two  different  grids  used  in  the  ZNSFLOW 
computations  produce  almost  identical  results  at  higher  angles  of  attack  and  are 
within  5%  of  each  other  at  lower  angles  of  attack.  The  computed  results 
predicted  by  the  CFD++  code  seem  to  match  the  data  slightly  better  except  at 
0  degree  angle  of  attack. 

Computatiojis  were  performed  for  the  H3P78  elliptical  projectile  at  a  higher 
Mach  number,  M  =  4,  and  several  angles  of  attack.  Figures  22  ttuough  25  show 
the  comparison  of  the  computed  aerodynamic  force  and  moment  coefficients 
with  the  data  at  M  =  4.  Figure  22  shows  the  comparison  of  normal  force 
coefficient  at  various  angles  of  attack.  The  normal  force  increases  with  the 
increasing  angles  of  attack,  as  seen  in  the  computed  and  measured  data.  The 
computed  data  were  obtained  with  Zeus  Euler  as  well  as  the  Navier-Stokes 
solvers  and  are  shown  in  open  symbols.  The  experimental  data  are  shown  in 
closed  circles.  All  computed  results  are  in  excellent  agreement  with  the  data  for 
all  angles  of  attack.  Computed  results  with  the  ZNSFLOW  code  again  were 
obtained  with  an  algebraic  Baldwin-Lomax  turbulence  model,  whereas  the 
CFD++  results  were  obtained  with  a  point-wise  two-ecjuation  turbulence  model. 
As  seen  in  Figure  23,  the  computed  pitching  moment  coefficients  predicted  by 
the  Euler  and  Navier-Stokes  techniques  are  in  excellent  agreement  with 
experimental  data.  A  comparison  of  the  computed  axial  force  coefficient  with  the 
data  is  shown  in  Figure  24.  As  expected,  the  largest  discrepancy  in  the 
comparison  is  with  the  Euler  code  prediction.  The  Euler  code  underpredicts  the 
axial  force,  as  seen  in  Figure  24.  The  axial  force  coefficient  predicted  by  the 
ZNSFLOW  code  matches  well  at  0  degree  angle  of  attack  and  is  within  10%  of 
the  measured  data  at  other  angles  of  attack.  Computed  results  obtained  with  die 
CFD++  code  are  generally  in  very  good  agreement  with  the  experimental  data. 
The  CFD++  code  was  also  used  to  compute  the  flow  fields  at  M  =  4  for  laminar 
flow  conditions. 


Figure  25  shows  the  comparison  of  laminar  and  turbulent  runs  wifli  the  data.  The 
laminar  results  at  all  angles  of  attack  are  xmderpredicted,  although  the 
discrepancy  seems  to  be  less  at  higher  angles  of  attack.  The  flow  in  the  wake 
region  behind  the  projectile  (and  tims  the  base  drag)  is  predicted  correctly  in  the 
turbulent  runs,  as  evidenced  by  the  very  good  comparison  with  flie  data. 
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Figure  18.  Computed  Particle  Traces,  M  =  4.0,  a  =  12. 
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Figure  19.  Normal  Force  Coefficient,  M  =  2.5. 


Figure  21.  Axial  Force  Coefficient,  M  =  2.5. 
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Figure  22.  Normal  Force  Coefficient,  M  =  4.0. 
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Figure  23.  Pitching  Moment  Coefficient,  M  =  4.0. 


5.3  Jet-On  Results 

Three-dimensional  calculations  at  Mach  4.0  and  several  angles  of  attack  between 
0°  and  12°  were  performed  for  the  same  eUiptical  projectile  (see  Figure  26)  for  the 
jet-on  conditions.  These  numerical  computations  were  performed  for  the  jet-on 
conditions  to  study  the  interaction  of  a  helium  jet  with  a  free  stream  M  =  4.0 
flow. 


Figure  26.  H3P78  Projectile,  Showing  the  Jet  Locations. 


The  real  gas  version  of  the  CFD++  code  was  used.  A  point-wise  two-equation 
turbulence  model  was  used  and  integrated  aU  the  way  to  the  wall.  Because  of 
symmetry,  only  one  of  the  two  jets  has  been  modeled  in  the  computations.  The 
ratio  of  the  jet  to  free  stream  pressure  was  varied  from  156  (100%  jet  pressure)  to 
39  (25%  jet  pressure).  The  velocity  of  the  jet  was  set  at  872  meters  per  second. 
Again,  the  projectile  geometry  and  a  set  of  experimental  wind  tunnel  data  for 
vahdation  of  the  computations  were  supplied  by  DERA,  UK  [23]. 


Figure  27.  Computed  Pressure  Contours,  Jet  On,  100%  Jet  Pressure,  M  =  4.0,  a  =  0°. 


The  computed  pressure  contours  in  the  symmetry  plane  are  shown  in  Figure  27. 
Here,  the  low  pressure  region  is  indicated  in  blue  and  high  pressures  are 
indicated  in  red.  The  red  region  on  the  lee  side  shows  the  effect  of  the  jet 
upstream  and  on  the  incoming  free  stream  flow.  As  expected,  the  pressure  in  the 
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base  region  or  in  the  near  wake  is  low.  The  flow  upstream  from  the  flare  junction 
remains  unaltered.  An  expanded  view  of  the  pressure  contours  in  the  base  region 
and  in  the  vicinity  of  the  jet  is  shown  in  Figure  28.  It  clearly  shows  the  effect  of 
the  jet  on  the  lee  side  flow  field.  The  pressure  field  on  the  wind  side  looks  similar 
to  that  of  the  jet-off  case.  A  circumferential  view  of  the  surface  pressures  on  the 
projectile  after-body  is  shown  in  Figure  29.  It  clearly  shows  the  strong  effect  of 
the  jet  in  all  directions  for  the  50%  jet  pressure  case.  The  surface  pressure  in  front 
of  the  jet  is  increased  while  the  surface  pressure  downstream  from  the  jet  is 
lower.  The  jet  affects  the  pressure  as  far  as  the  centerline  of  symmetry,  which 
indicates  the  strong  interaction  of  the  two  jets  in  that  area.  A  longitudinal  cut  at 
the  jet  centerline  was  selected  and  used  to  show  the  concentration  of  helium  in 
that  plane  containing  the  jet  centerline  (see  Figure  30)  for  different  jet  pressures. 
Here,  black  represents  zero  helium  concentration  (100%  air)  and  yellow 
represents  high  helium  concentration.  At  the  jet  exit,  flie  helium  concentration  is 
one.  This  figure  clearly  shows  the  helium  jet  interacting  with  the  free  stream 
computed  surface  pressure  contours  for  the  three  different  jet  pressure  cases. 
Computed  surface  pressures  have  been  obtained  along  the  centerline  of 
symmetry  as  well  as  along  the  centerline  of  the  jet  itself.  These  surface  pressures 
are  used  to  study  the  effect  of  the  aerodynamic  interference  resulting  from  the  jet 
interaction  with  the  free  stream  flow.  Computed  surface  pressures  have  been 
compared  with  the  experimental  data.  Figure  31  shows  the  comparison  of  the 
computed  surface  pressures  with  the  data  measured  along  the  centerline  of 
symmetry  for  various  jet  pressure  cases.  The  computed  results  here  are  shown  as 
symbols  and  the  data  are  shown  as  lines. 


Figure  28.  Computed  Pressure  Contours  Near  the  Jet  and  Wake,  Jet  On,  100%  Jet 
Pressure,  M  =  4.0,  a  =  0°. 
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Figure  29.  Computed  Surface  Pressure  Contours  Near  the  Jet,  50%  Jet  Pressure, 
M  =  4.0,  a  =  0°. 


CFD  and  the  data  show  a  pressure  rise;  however,  the  numerical  predictions 
indicate  a  much  stronger  jet-to-jet  interaction  than  is  indicated  by  tiie  data.  The 
pressure  rise  is  somewhat  decreased  with  decreasing  jet  pressures.  The 
computed  results  presented  here  used  the  fine  mesh  (grid  3),  and  it  is  believed 
that  the  discrepancy  is  not  attributable  to  the  number  of  mesh  points  used  in  tiie 
computations. 

The  comparison  of  the  computed  surface  pressures  along  the  centerline  of  the  jet 
for  various  jet  pressure  cases  is  shown  in  Figure  32.  Again,  computed  results 
have  been  compared  with  the  available  data.  Except  for  the  100%  jet  pressure 
case,  the  computed  results  are  in  good  agreement  with  tire  measured  data.  CFD 
and  the  data  clearly  show  the  expected  pressure  rise  ahead  of  the  jet.  Also,  this 
figure  shows  the  lower  pressures  downstream  from  the  jet  (pressure  less  than  tire 
free  stream  static  pressure).  Although  not  shown  here,  this  level  of  lower 
pressure  downstream  from  the  jet  is  also  very  similar  to  that  measured  in  the 
experiment.  All  tiie  flow  field  effects  associated  the  jet  interaction  have  not  yet 
fuUy  been  analyzed.  Overall,  the  computed  flow  field  around  the  projectile  and 
the  jet  showed  similar  characteristics  as  those  observed  in  the  experiment. 
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(a) 


(b) 


(c) 


Figure  30.  Computed  Helium  Concentration  Contours  Near  the  Base  Region, 
(a)  100%,  (b)  50%,  and  (c)  25%  Jet  Pressures,  M  =  4.0,  a  =  0°. 
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re  31.  Surface  Pressure  Contours  Comparison  Along  the  Centerline  of 
Symmetry,  M  =  4.0,  a  =  0°. 
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Figure  32.  Surface  Pressure  Contours  Comparison  Along  the  Centerline  of  Jet, 
M  =  4.0,  a  =  0°. 


6.  Conclusion 


CFD  approaches  were  used  to  compute  the  supersonic  and  hypersonic  flow 
fields  and  aerodynamic  forces  and  moments  on  elliptical  projectiles.  Steady  state 
numerical  results  have  been  obtained  at  several  supersonic  Mach  numbers 
between  2.5  and  4.0  and  several  angles  of  attack  from  0°  to  12°  for  die  jet-off 
conditions  with  the  use  of  Euler  and  Navier-Stokes  flow  solvers.  In  addition, 
numerical  computations  have  been  performed  for  the  jet-on  coriditions  to  study 
the  interaction  of  a  helium  jet  with  a  free  stream  M  =  4.0  flow.  In  general,  very 
good  agreement  of  the  computed  aerodynamic  coefficients  with  the  experimental 
data  was  achieved  at  all  speeds  and  angles  of  attack  investigated  for  jet-off 
conditions.  Both  Navier-Stokes  codes  predicted  the  normal  force  and  pitching 
moment  coefficients  very  well.  However,  the  CFD++  code  seems  to  predict  the 
axial  force  more  accurately,  partly  because  of  the  advanced  turbulence  modeling 
used  in  the  computations.  CFD  results  for  the  jet-on  case  show  that  the 
qualitative  features  and  strong  flow  interaction  between  the  jet  and  die  free 
stream  flow  are  similar  to  those  observed  in  the  experiment.  Computed  surface 
pressures  along  the  jet  centerline  compare  better  than  diose  along  the  line  of 
symmetry.  The  results  show  the  predictive  capabilities  of  CFD  techniques  for 
supersonic  flow  over  elliptical  projectiles  with  and  without  jet  interaction. 
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